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AbstractE 

Imagine that measurements are made at times to and t\ of the trajectory of a physical 
system whose governing laws are given approximately by a class A of so-called prior vector 
fields. Because the physical laws are not known precisely, it might be that the measurements 
are not realised by the integral curve of any prior field. We want to estimate the behaviour 
of the physical system between times to and t\. 

An integral curve of an arbitrary vector field X is said to be feasible when it interpolates 
the measurements. When X is critical for 1? distance to A, the feasible curve is called 
a conditional extremum. When the distance to A is actually minimal, the conditional 
extremum is a best estimate for the intermediate behaviour of the system. 

The present paper does some of basic groundwork for investigating mathematical prop- 
erties of conditional extrema, focusing on cases where A comprises a single prior field A. 
When A = a conditional extremal is a geodesic arc, but this special case is not very rep- 
resentative. In general, A enters into the Euler-Lagrange equation for conditional extrema, 
and more can be said when A is conservative or has special symmetry. 

We characterise conservative priors on simply-connected Riemannian manifolds in terms of 
their conditional extrema: when time is reversed, a constant is added to the L 2 distance. 
For some symmetric priors on space forms we obtain conditional extrema in terms of the 
Weierstrass elliptic function. For left-invariant priors on bi-invariant Lie groups, conditional 
extrema are shown to be right translations of pointwise-products of 1-parameter subgroups. 

^010 MSC: Primary: 49K15, 53C22, 37M99; Secondary: 37C10, 37E35 
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1 Introduction 



This paper focuses on the following question. 

Let A be a given C°° vector field on a connected C°° m-manifold M . A curve x in M is 
thought to be an integral curve of a unknown vector field that is near A. Precise observations 
x(tk) are made of x at finitely many times t^ where to < t% < . . . < t n , but no further 
information about x is known. How then to estimate x(t) for t ^ t^? 

Consider first the case where M is Euclidean m-space E m . If n = 1 and t £ (to,ti) we 
could use the weighted average 



x(t) 



(tl - t)XQ + {t- t )xi 

ti — to 



or, more generally, for n > 1 the natural cubic spline [2]. But these interpolants take no 
account of the information that x^ ~ A(x), where x^ denotes the derivative of x. To do 
so requires knowledge about the nature of the supposed approximation, which in practice 
depends on context. We are going to assume x^ 1 ' is close to A[x) in the I? sense, namely 
x W(t) = A{x(t)) + g(t) where g : [t ,t n ] -)• E m is L 2 -small. Possibly g(t) = G(x(t)) where 
G is another vector field. 

In practice there might be a class A of vector fields rather than a single prior field A, 
with the information that x^ ~ A(x) for some A G A. When A is large the quality of 
the information about the field is small. For A the space of all C°° vector fields there is 
no information except the observations, and we are reduced to classical methods of inter- 
polation. Between these extremes, where A is a singleton or very large, A might instead 
be parameterised by a finite-dimensional a G W. The literature on model building and 
parameter selection for dynamical systems (see [B], pQ) includes cases where observations 
are contaminated by noise, as in [8]. In the simpler setting of our question, there is no 
noise, and the requirement is for interpolation rather than approximation. An important 
case for the present paper is where A is a singleton, and parameter selection is not an issue. 

Our task is to study nonparametric interpolants x minimising the 1? difference between 
x^ and A[x) subject only to x{tk) = Xk for < k < n - a problem in the calculus 
of variations. The present work also finds new results and explicit solutions for systems 
with symmetry in the setting of Riemannian geometry, where E m is replaced by an Tri- 
dimensional Riemannian manifold M. 



2 Description of Results 

In ^3] conditional minima and extrema x are defined relative to a class A of vector fields. 
Working first within the class of so-called almost-C 2 curves, it follows from Theorem[l]in 
that if A is C k then, except at t\,t2, ■ ■ ■ ,t n —\, a conditional extremum x is automatically 
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C k+l . Examples [Xj § g point out some simple facts. For instance, when .A = {0} a 
conditional extremum is the same as a geodesic: conditional extrema relative to prior 
fields generalise geodesies. 

Theorem [I] also gives the Euler-Lagrange equation ^ for a conditional extremum x with 
respect to a prior field A. The left hand side of ^ is the covariant acceleration of x. The 
right hand side has two nontrivial terms, namely the gradient of the squared norm of A, 
and a second term made by 

• replacing A by the exterior 1-form A T corresponding under the Riemannian metric 

• taking the exterior derivative of A T 

• contracting with x^ 1 ' 

• replacing the resulting 1-form by the corresponding vector field defined along x. 

A conditional extremum is a C° track -sum of solutions of Q . Consequently, when A is a 
singleton it suffices to consider the case n = 1. For M = E m , Example [4] solves equation 
@ when A is an affine vector field (Example [3] finds x directly when A is constant). 
Theorem [2] shows, for an arbitrary Riemannian manifold M, that 1 1 a?^ 1 ) ) 1 1 2 — ||t4(x(£))|| 2 
is conserved along solutions x of Q, and that if A is bounded and M is complete then 
x extends to a solution of Q defined on the whole of R. This generalises the well-known 
results for geodesies, that ||xW|| 2 is conserved, and that complete Riemannian manifolds 
are geodesically complete [5] Part II. 

For a geodesic x : [io 5 *i] — > M its reverse x : [— to] — > M, given by x(u) = x(—u) is 
also a geodesic. From the asymmetric role of the prior field A in the definition, it seems 
improbable that the reverse of a conditional extremum would be a conditional extremum 
for the time-reversed data. The improbable happens when the 1-form A T is closed, as 
proved in Corollary [T] of Theorem [TJ The reverse of a conditional optimum need not be 
optimal, even when A T is closed, as seen in Example^ 

A stronger condition is for A to be reflexive, namely J(x, A) differs from J(x, A) by a con- 
stant depending only on A and independent of x. Theorem [3] of ^5] shows A is reflexive if 
A T is exact and, conversely, if A is reflexive then A T is closed. So if M is simply connected 
then reflexive is equivalent to conservative, as illustrated in Example [7] by numerical cal- 
culation. These kinds of calculations are performed by replacing curves x by lists of points 
in W where p is 2, 3 or 4, then numerically optimising with Mathematica's FindMinimum. 
In applications one might rely on refinements of such methods, but the present paper gives 
theoretical results, including closed form solutions for conditional optima, when M and A 
exhibit symmetry. 

In M is a semisimple Lie group G with bi-invariant Riemannian metric. Whereas 
geodesies in G are right-translations of 1-parameter subgroups, Theorem [4] says solutions 
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of ([2]) for a left-invariant prior field A are right-translations of pointwise products of a pair 
of 1-parameter subgroups. One subgroup is generated by the value of A at the identity 
of G. The other is obtained from the first after comparison with x$ and x\. Corollary 
[4] extends this to n > 1. So the Euler-Lagrange equation ^ is solved in closed form, as 



illustrated for the group S 3 = SU(2) of unit quaternions by Example 10 All this supposes 
that the left-invariant field A is given, as it would be if A was a singleton. For a finite set 
of left-invariant fields, J is minimised for each A, and the smallest is chosen. 

For A an infinite set of left-invariant vector fields on G, Corollary [4] of Theorem [4] reduces the 
search for an optimal pair to a finite-dimensional optimisation problem, whose numerical 
solution is straightforward. We also prove some theoretical results for the case when A is 
generated by a submanifold of Q. Theorem [5] gives a nonlinear equation for A in terms 
of the exponential map of G. When A is the set of all left-invariant fields and the tk are 
equally spaced, the equation takes a simpler form given in Corollary |6j and verified in 
Example 11 Specialising to n = 2 (three observations), Corollary [7] gives a simple solution 
for (x,A) in terms of the exponential. 

From ^7| onwards, M is either the unit sphere S 2 in E 3 or the unit two-sheeted hyperboloid 
H 2 in Lorentz 3-space, the two cases being treated simultaneously. The prior fields A, 
parameterised by a pair of functions ft, 7 : R — >• R, are invariant with respect to rotations 
in M 2 x {0}. Under these conditions, in equation ([2]) is rewritten as three coupled 
nonlinear 2nd order scalar ODEs Q, ([7]), ([8| for the scalar- valued coordinates x\, X2, X3 
of a conditional extremum x. As well as equation ^ which reappears as Q, rotational 
symmetry gives another conserved quantity (10). Then x is found by quadrature in terms 
of X3 : R — > R, and a first order ODE (|12j) for X3 is given. Then ft and 7 are taken as 
constant, and the solutions of ([7]), Q for which X3 is constant are found. Further 
results for X3 nonconstant depend on whether A has a longitudinal component, namely 
whether 7 7^ 0. 

For 7 = 0, Theorem [6] of ^solves 0, ^ for x in terms of sin, cos, sinh, arctan, cosh 
and arctanh. For both S 2 and H 2 there are two different kinds of solution, depending on 



whether (0,0, 1) lies in the image of x. Examples 13, 14 describe solutions where M is S 2 
and H 2 respectively, with H 2 replaced by the Poincare unit disc for ease of illustration. 
As shown in Figure |4j there exist non-optimal conditional extrema. 

In §7.2[ 7^0 and Theorem [7] gives £3 in terms of the Weierstrass elliptic function p. In 
particular X3 is periodic, and solutions x of Q, 0, Q are found by quadrature in terms 
of p. Even when A is conservative, namely ft = 0, these solutions can be geometrically 
interesting as shown for M = S 2 in Figure [8j When A is not conservative A is usually 
not reflexive (never on S 2 or H 2 ), and Example 16 describes a conditional extremum 
x : [0, 1] — > H 2 and its reverse x where A is nonconservative. This is illustrated in Figure 
[9] by mapping into the Poincare disc. 
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3 Optimality 



Let ( , ) be a C°° Riemannian metric on M. For ao < a\ G E, a continuous curve 
y : [aO) a i] — > M is said to be almost-C 2 when, for some a > 0, some ao = so < *i < 
. . . < s p = ai and all j = 1, 2, . . .p, each restriction y| [sj-i, Sj] extends to a C 2 curve in M 
defined on (sj-i — 5, Sj + 5). Then so,s±,... s p are singular parameters of y. 

For some n > 1, let xo, xi, . . . x n G M and to < h < • • • < t n G M be given. An almost-C 2 
curve x : [to, t n ] — > M is said to be feasible when x(tk) = Xk for all k = 1, 2, . . . , n. The set 
of all feasible curves is denoted by X = Xt ,t 1 ,...,t n ;x ,x 1 ,...x n - Given also a nonempty set A 
of C 1 vector fields Aon M, define J : x A ->• E by 

J(rr,A) := / ||xW(t) - A(x(t))|| 2 dt 

where x^ 1 ^ is the derivative of x with respect to t, and || || denotes the Riemannian norm. 
A minimizer (x, A) of J = Jt ,t 1 ,...,t n ;x ,x 1 ,...x n is said to be conditionally- optimal or just 
optimal. When there is no doubt about A, for instance when A is a singleton, we say x is 
conditionally optimal rather than (x,A). 

Suppose there is an unknown curve z : [to, in] - ► approximately satisfying 

z«(t) = 

for some vector field A in a given parameterized set A of prior fields. To say (x, A) is 
conditionally optimal means the velocity field of the feasible curve x is as near as possible to 
A £ A while agreeing with observations made from z at the parameter values to, t\, . . . , t n . 

Example 1 If some x G X is an integral curve for some A £ A then (x, A) is conditionally 
optimal. For n = 1 and xo = x\, the constant curve is not necessarily optimal. □ 

Example 2 Let A = {0}. Then x G X is conditionally optimal when it is a track sum of 
n minimal geodesic arcs from the Xk-i to x&. In particular, when n = 1, x is a minimal 
geodesic arc. □ 

Example 3 Let M be Euclidean m-space E m , and let A be a constant field. Then 

J{x,A) = f n \\ x ( 1 )(t)\\ 2 dt-2(A,x n -x } + {t n -t )\\A\\ 2 . 
J t 

So, for any nonempty set A of constant vector fields on E m , a necessary condition for (x, A) 
to be conditionally optimal is that x be the piecewise-affine interpolant of the observed points 
Xk at the tk- As will be seen, it is rare for x to be independent of A in this way. 
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In the least informative situation, where A is the set of all constant fields there is a unique 
optimal pair (x,A), with A = (x n — xq) / {t n — to) ■ It might seem strange that the optimal A 
takes no account of x±, X2, ■ ■ ■ , %n—l> but the observations are noise-free. The uncertainty 
concerns only x(t) for t 7^ So A should indeed be a weighted mean of intermediate 
estimates: 

*-i>-^>(^r)- 

The outcome is less predictable in Example ^ where A is the set of all affine vector fields 
on E m , and in {6.1 where A is all left-invariant vector fields on a semisimple Lie group. 
□ 



4 First Order Necessary Conditions 

Let x : [to,*n] x ( — e i e ) ~~ > M be continuous where e > 0. Set Xh(t) := x'(/i) := x(t,h). 
Given x G X with singular parameters sq, s\, . . . , s p , we call x a variation of x when 

• x(t, 0) = x(t) for all t G to, t n ] 

• Xh G X for all h G (— e, e) 

• for some a > 0, each restriction x\[sj-i,Sj] x (— e, e) extends to a C 2 map from 
(sj-i — a, Sj + a) x (— e, e) to M. 

So each Xh has the same singular parameters as x, and x t is C 2 for all t G [toj^n]- Set 
W(t) = (x*)'(0) where ' denotes differentiation with respect to h and, for s G (to,t n ), 
Ax^(s) := lim x^'{u) — lim a;W(u). Then Ax^(s) = for s nonsingular. A vector 

field X on M corresponds to a differential 1-form given by X T (Y) := (X,Y), where 
Y is any vector field. Of course A is conservative precisely when A T is exact. Since A G A 
C 1 so is A T , and another differential 1-form #a,x on M is given by 

^ -> 0A,xPO = -0A,yPO := (cL4 T )(X,Y) = 

I((i,F))-F((AI})-(A,[I,F]) = (VxA^)-(VyA^) 

because the Levi-Civita covariant derivative V is torsion-free. The differential 1-form 9a,x 
corresponds to another vector field Q~ A T X - From the definition of 9a,x, &a,x(X) = —6a.y{X) 
and, for any C 1 function a : M — > M, 

W = a^A,x + X(a)^ T -(A^)da. (1) 
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Lemma 1 For any variation x of x E X , 

l -^J( Xh ,A)\ h= , = P jyW{ Sj ),^\ Sj ))- f\w,V t x^ dt. 



Proof: Because V is symmetric and compatible with the Riemannian metric, 

±JLj(x h ,A) = f\v t x' h -V h A,x^\t)-A(x h {t)))dt 
where ' denotes differentiation with respect to h. On integration by parts this becomes 



X>U*;), A4 1} ( Si )) - f\xlV t {x^\t) - A(x h (t)))) + (V h A,x£\t) - A(x h (t))) dt. 



3=1 



to 



1 d 

Setting /i = 0, - — J(x/ l , 71)1^=0 = 

P ^(W( S] ),Ax^( S] ))- r n (^(t),V t (x( 1 )(t)-A(x(t)))) + (V^(x(t)),x( 1 )(t)-^(x(t)))(it 

3=1 Jt ° 

p-l 

'W{ Sj ),Ax^\s ))- 

3=1 Jto 
/-I 

[W( Sj ),AxU( Sj )}- 

3=1 Jto 
□ 

We say (x, A) E X x A is J -critical (or just critical) when, for all variations x of x, 

—7(^,^4)1^=0 = 0. 
Then x is called a conditional extremum. Prom Lemma [T] follows 

Theorem 1 If{x, A) E XxA is critical if and only ifx has singular parameters to, ti, . . . i r 
and, /or all k = 1,2, ... ,n and all t E (ifc-i, ifc), 



Ax^( Sj )) - f\w,V t x^(t) - ^gradPII 2 )) - 9 AxW (W) dt 
J2(W( Sj ),AxW( Sj )) - j\w,V t x^(t) - ^gradP|| 2 - e£ w ) dt. 



V t x«(i) = ^gradP|| 2 + ^ (1) . (2) 



□ 



In order to be optimal, J should also be critical with respect to variations in A. In the 
present paper we assume that A has been found somehow; perhaps A is a singleton. By 
Theorem [TJ a conditional extremum is a C 2 track -sum of feasible curves satisfying »2v, 
namely conditional extrema for the case n = 1. 

If A = 0, a conditional extremum is the same as a geodesic. 

Example 4 Let n = 1 and M = E" 71 . Identifying vector fields on E m with functions 
A : E m — > £ m ; 

M = (^A(x(t)) J y)-( < M.(y),x«>. 

Ze£ ^4 6e affine, of the form A(y) = By + c where c G and -B is linear. Then x is 
critical when, for all t G (to,ti), 

X ( 2 )(t) = B T (Bx(t) + c)+SxW(t)-B T xW(t) 

where T means matrix transpose. Equivalently z^\t) + S T z = S T c where z(t) : = 
x^(t)-Bx{t). So z(t) = e- wT d + c where deE m , and 

x (t) = e ^ B yo + e tB f e~ sB (e~ sBT d + c) ds 

where d is chosen so that x(t\) = x\. Example^ dealt with B = 0. □ 

Theorem 2 Let x : [to,ti] M be a solution of For some b G K and a// f G (ioj^l); 

„(i)mi|2 = ||A(a:(t))|| 2 + 6. (3) 



// M zs complete as a metric space and A is C°° and uniformly bounded, then x extends 
to a unique C°° solution of defined on all ofM.. 

Proof: Taking inner products with x^ of both sides of 

2dt {X ' X > ~ {VtX ' X } " 2di" 11 + ^>^ x ) ~ 2dt m ' 

Integrating both sides, Q follows. Alternatively, Q follows from Noether's theorem and 
time-invariance of the Lagrangian. 

Suppose now that M is complete, and that A is C°° and uniformly bounded. For some 
real /3 and all y G M we have ||^4(y)|| 2 < f3 where || || is the Riemannian norm. By the 
Picard theorem on solvability [4 J §7.5, x is C°°. By m, \\x^\\ 2 <f3 + b. These statements 
hold also for any C 2 extension of x to a solution of ( 2 ) . 
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Let 7+ (respectively 71) be the nonempty sets of real t\ > t\ (respectively to < to) such 
that x extends to a C°° solution of ^ defined on (to, ii) (respectively (to, t±). To complete 
the proof it suffices to show that 7+ is not bounded above and 71 is not bounded below. 

If 7+ is bounded above, let t\ = sup7+. Then, for any integer p > P with P sufficiently 
large, x extends to a C°° solution x^ : (to,ti — l/p) — > M of Q. Set y p := x^(i\ — 2/p). 

Because the II^njII are uniformly bounded for p > P, {y p : p > P} C M is Cauchy, with 

limit y say. Using some coordinate chart containing y, represent each x^(ti — 2/p) as a 

vector z p in E m . Because the i?) are bounded with respect to the Riemannian norm, the 
sequence {z p : p > P} C E m is also bounded in the Euclidean norm. So {z p : p > P} has 
a convergent subsequence {z Pq : q > 1} whose limit represents z G TMg. 

By the Picard theorem, for some r > 0, any i 6 R, and any (y, 5) G TM sufficiently near 
(y, z), there is a unique C°° solution x : (i — r, £ + r) — >• M of d2| satisfying x(£) = y and 

xW(t) = i. 

Choose q so large that, with t = t\ — 2/p g , 

• (y,^) := ( ;E [Pg](*)' a; [ J !, ) ](^)) is sufficiently near (y, z) and 

• 2/p q < T. 

Splicing X[ p ] and x gives a C°° extension of x defined over (to,ii + T — 2/p q ). Since 2/p q < r 
this contradicts the definition of t\ as sup 7+ . The proof that 71 is not bounded below is 
entirely similar. □ 

Example 5 For ^4 = the hypothesis that M be complete is needed for extendability, 
from the Hopf-Rinow Theorem. Extendability may also fail when M is complete and A is 
unbounded, as when M = E 1 and A(x) = x. □ 

5 Closed Prior Fields 

The definition of J(x, A) depends not only on x G X and A G A, but also on the t(.£l 
and the xt G M which are used to define X . The reverse x of x is defined by x(u) := x(—u) 
for u G [— t n , —to] and 

X := X-t n ,-t n -i,...,-to;x n ,x n -i,...xoi J := J-t n -t n -i,...,-to;x n ,x n -i,...xo ■ X X A R. 

If x G X then From Theorem [l] follows 

Corollary 1 Suppose dA T = 0. Then (x, A) G X x „4 is J -critical if and only if (x, A) G 
<Y x „4 is J -critical. 
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Proof: dA T = ==> 9 A x (i) = = 6 A , and so (x, A) is critical if and only if 
V t 3«(t) = ^^(Pll 2 ) ^ V t -i«(t) = ^grad(||A|| 2 ). 

□ 

However if (M T = and (x, A) is J-optimal, (x, -A) need not be J-optimal. 

Example 6 Let A be the clockwise unit vector field on M = S 1 with the standard Rie- 
mannian metric. Take n = 1, to = 0, t\ = 2n and xq = x\ = (1,0). Then (x,A) given by 
x(t) = (cost, — sin t) is J-optimal, with J(x,A) = 0, and J(x,A) = 8ir. Also (x*,A) given 
by x*{u) = (cosm, — sinu) for u 6 [— 2ir, 0] is J-optimal, with J(x*,A) = 0. □ 

Definition 1 A nonempty set A of C 1 vector fields on M is said to be reflexive when, 
for any real to < t± < . . . < t n , and any xo, x%, . . . , x n £ M, there exists CA, Xo ,x n G K 
independent of to, t%, . . . , t n and independent of xi, X2, ■ ■ ■ , x n -i, and there exists A £ A 
depending only on A, such that, for all x £ X , J(x,A) = J(x,A) + ca- Call A reflexive 
when {A} is reflexive. □ 

In particular, if A is reflexive then for all (x, A) £ X x A 

(x, A) £ X x A is J-critical <^=^ (x, A) £ X x A is J-critical 

and 

(x,A) £ X x A is J-optimal <^=^ (x, A) £ X x A is J-optimal. 

If, as in Example [2j A is closed under multiplication by —1 then A is reflexive, with 
c A,x ,x n = and A = —A. In Example [3] any nonempty set of constant vector fields on E m 
is reflexive, with ca,xo,x„ 

, x n — xq) and A — A. The condition that A be reflexive is 
especially stringent when A is a singleton. 

Example 7 For M the unit 2-sphere S 2 in E 5 , set A = {A} where A(v\,V2,V3) = 
(v2,— Vi,Q) for v £ S 2 . The integral curves of A are shown (black) as latitudinal cir- 
cles traversed in the clockwise direction. Set n = 1, to = 0, t\ = 1, xo = (0.866,0,0.5) 
and Xl = (0.5187,0.8486,0.1039). Figure [7J sh ows (red) the conditional minimum x and 
(blue) the conditional minimum x* for the reverse data, as approximated from a numer- 
ical computation. Theorem [6] gives closed form expressions for such curves. We find 
J(x*,A) ps 0.18 < J(x,A) « 0.44. So A is not reflexive. □ 

Indeed in Example [7] A cannot be reflexive, by the first part of 

Theorem 3 If A is reflexive then dA T = 0. // A T = d(f> where 4> : M — > K then A is 
reflexive, with CA, X0 ,x n = 4((f>(x n ) - 4>(x )). 
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Figure 1: Conditional Minima With A Not Reflexive, in Example [7] 

Proof: Suppose first that A is reflexive, and choose to < t\ < . . . < t n and xo, xx,...,x n £ 
M. For any x G X, and any variation x of x, 2K(h) := J(xh, A) — J(x~h, A) = ca,x ,x u - By 
Lemma [TJ = K' (0) = 



p-i 



j'=i 



^(^( Sj ),Ax«( Sj )>- / {W,V t x^{t)--g^d\\Af-e 



p-i 



3=1 



to 



-to 



dt 



^(Wi-s^Ax^i-sj)) + / (W(t),VtxU(t) - -gradpf -9 



A,xW' 



dt 



to 



... 

(W(t),V t x^(t)--gr a d\\A\\ 2 -el x(1) )dt+ I {w{t),v t x^\^--g^d\\A\\ 2 +e- T xW )dt 



to 



2 1 (W(t),9- T xW ) dL 
to 



Since this holds for all variations x, we have &a xw\x(t) = f° r an nonsingular t and, by 
continuity, for all t E [to, in]- Since this holds for all x £ X, dA T = 0. 



Suppose now that uja = d(f>. Then J(x, A) — J(x, A) 



to 



ftn 

\x^(t) - A(x(t))\\ 2 - \\x^(t) + A(x(t))\\ 2 dt = -4 / {A(x(t)),xW(t)) 

Jto 



dt 



{<t>ox{t)) dt = - 4(cf>(x n ) - <f>(xo)). 



to 



□ 
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Corollary 2 Let H 1 (M;R) = 0. Then A is reflexive if and only if A is conservative. □ 

Example 8 For p > 2 let M be any Riemannian manifold with the homotopy type of 
the unit sphere S p in E p+ , or real projective p-space RP P , or complex projective m-space 
CP P . Then A is reflexive if and only if A is conservative. In particular, regardless of the 
Riemannian metric on M = SO(3), A is reflexive if and only if A is conservative. □ 

6 Left-Invariant Priors on Bi-Invariant Lie Groups 

Take M to be a Lie group G with a bi- invariant Riemannian metric ( , ), and let A be a 
left-invariant vector field on G [5], [3], |10| . Let n = 1. The left Lie-reduction of a vector 
field Z defined along x : [t$,ti\ — > G is defined to be the curve Zl in the Lie algebra 
Q = TG\ given by 

Z(x(t)) =: dL(x(t))i(Z L (t)) 

where L is left-multiplication, and 1 is the identity element of G. Then Al is constant. 
Denote (x^)l by V. Because ( , ) and A are left-invariant, for any Y G TG x a\, 

(Yl,(0- T xW )l) = (Y,e- T x(1) ) = 9 AxW (Y) = (dA T )(x^,Y) = (V t A,Y)-(V Y A,xW) = 

((V t A) L ,Y L )-{(V Y A) L ,V) = ±([V,A L ],Y L )-±([Y L ,A L ],V) = {[V,A L ],Y L ) 

namely = [V, Al]- Take n = 1. Then a feasible curve x is a conditional extremal 

when V (1) (t) = [V(t),A L ]. Equivalently V(t) = Ad( e - (i -' (,)Ai )F(t ), and 

dRie-^^^odLix^odLie-^-^ 

= dL(x(t)e~( t ~ t °^ AL )i(V(to) — Al). 

So x{t)e-^- to)AL = e (*-*o)(^(*o)-A z ,) < Taking B L := V(t ) - A L , this proves the first 
assertion of 

Theorem 4 Let A be a left-invariant vector field on a Lie group G with a bi-invariant 
Riemannian metric. Let n = 1. A feasible curve x : [tci^i] G is a conditional extremum 
if and only if, for all t G [to, ti], 

x (t) = e {t - to)BL e {t - to)AL x (4) 

where P>l G Q is such that x{t\) = x\. Then, for all t G (to,ti), 

\\x^(t)\\ = \\A l +Bl\\, \\x (1 Ht)-A(x(t))\\ = \\B L \\, (xW(t),A(x(t))) = (A L +B L ,A L ). 
In particular J(x,A) = (t\ — to)||^i|| 2 - 
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Proof: Let x be a conditional extremum. Because ( , ) is left-invariant, Hx^^i)]! = 
||V(t)|| = ||V(to)||, since Ad acts by isometries. So ||a;W(t)|| = \\A L + B L \\. Because A is 
left-invariant and V(t) = Ad(e'^ to ^ AL )V(t ), 

\\V(t)-A L \\ = \\Ad(e-^ AL )V(t )-A L \\ = \\V(t ) - A L \\ = ||B L || 

because the isometry Ad(e^~ to ^ AL ) : Q — > Q fixes Al- Similarly, (x^(t),A(x(t))) is 

(V(t),A L ) = (Ad(e-( t - t ^)V(t ),A L )= (V(t ),A L ) = (A L + B L ,A L ). 

□ 

Example 9 Suppose x\ = e t,AL xo for some t* > to- Substituting Bi = — -Al in 

ti — to 

we find that the geodesic given by 

(,t-t o xu-t ) . 
x(t) = e 'i-*o xo 

is a conditional extremum. □ 

Corollary 3 Let A be a left-invariant vector field on a bi-invariant Lie group G. Then 
x : K — > G satisfies equation |§|) with x(0) = 1 if and only if x is a pointwise product 

x (t) = e tBL e tAL 

of one-parameter subgroups of G. □ 

Example 10 Take G to be the group S 3 C H of unit quaternions, with bi-invariant Rie- 
mannian metric from E 4 = EL A left-invariant vector field A corresponds to a pure imagi- 
nary tt£H, where A(x) = xa. The one-parameter subgroups of S 3 are given by gp(t) = e t ^ 
where f3 G H is pure imaginary. Choose q a ,q/3 6 S 3 so that q a aq a = ai and qpfiqp = bi 
where a,b 6R and q is the conjugate of q £E. Solutions x : K — > S 3 of with x(0) = 1 
have the form 

x(t) = q a (cos(at) + i sm(at))q a q~p(cos(bt) + i sin(bt))qp. 
Since A is nowhere-zero and M is compact, A is not conservative. 

TakeA L = (-0.5, -0.5, 0.3), x = 1 andx x = (-0.0359448,-0.228089,-0.937324, -0.260972). 
Figure^ shows (red) exp" 1 ox: [0, ir] — > E 3 for a conditional minimum x : [0, 7r] — > S 3 . 
The preimage under exp of the orbit y : [0, n] — > S 3 of A is shown (black) as a directed line 
segment, with (blue) the preimage oft \-> e tBL y(it) fort 6 [0, tt], where Bl = (0.2,0.2,0.2). 
□ 

For 1 < k < n set s& := t& — t^-i and y& := XkxZ_y 

13 
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Figure 2: A Pointwise Product of 1-Parameter Subgroups of S" 3 , in Example 10 



Corollary 4 Let (x,A) be optimal, where A is a left-invariant vector field on a Lie group 
G with bi-invariant Riemannian metric. For some Bli, Bl,2, ■ ■ ■ , Bi, n G G> 

x (t) = e (*-**-i)^.» e (*-t*-xM ia . fc _ 1 for t€[t k -!,t k ] 
where e SkBh ' k e SkAh = y\~, and the geodesic arcs s G [0, Sk] \— > e sBl " k G G are minimal. □ 

When .A is finite, Corollary [4] allows us to compute all optimal pairs (x, A) in terms of the 
exponential map of G. In particular, if for each k there is a unique minimal geodesic on 
[0, Sk] joining e SkAh to y^, then x is uniquely defined by A. 

6.1 Manifolds of Left-Invariant Priors 

Let A be a set of left-invariant vector fields on G where Al '■= {A^ : A G A} is a C 1 
immersed submanifold of Q of dimension p < m. The following analysis can be developed 
for when Al is the image of a C 1 immersion (not necessarily one-to-one), but we are 
interested in cases where Al is an affine subspace of Q, especially Al = Q- 

Theorem 5 Let (x,A) G X x A be optimal, where e~ SkAL yk is not a critical value of the 
exponential map for any k. Defining Bl^ i n terms of Al, as in Corollary^ 



l _ e -s k ndA L 



L.k 



k=l 

is orthogonal to the tangent space of Al at Al. 
Proof: For Al near Al, the relation 

e s k B L , k _ yk e -s k A L ^ 
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determines Blu locally as a C°° function of Al, with value Blu at Al- Using the formula 
for the derivative of the exponential map [9], [TO] to differentiate ^ at Al, 

^(e'*^.*) o f J o dB Ltk (W) = dL(y k ) o dL( e - s ^) o f — j (W) 

where L(g) is left multiplication by g G G, and W is tangent to .At, at A^,. Then by ([5]), 

1 _p-SfeadB ijfc \ /-■ _ s k adA L \ 

(B £)fc ,dB L , fc (W0) = (Si, fc ,f — Jodfl £ , fc (W)) = 

/ 1 _ pSfeadA^ \ /-i _ -s k adA L \ 

By Theorem|4j Al is a critical point of J(x,A) = J2k=i s k\\^L,k(^L)\\ 2 ■ So 
for all tangent to Al at Al. □ 



Corollary 5 Let A L = Q. Then (Al, V] s k B L ,k) = 0. 

fc=i 

Proof: Since B L = 0, = (A L , Sfc ^ J 5 L 

"2^^^ — J A L,B L<k ) = -2^s fc (AL 5J B L) 



□ 



Corollary 6 Let Al = Q, s\ = S2 = ■ ■ ■ = s n and suppose s\Al is not a critical point of 
the exponential map. Then, with the hypotheses of Theorem^ 

n 

J2 B L,k = o. 

fe=l 
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Proof: The derivative of the exponential map at s\Al is dL(e SlJ 



si ad 



□ 



Example 11 Let G = S 3 , Al = Q , n = 4 and s% = S2 = S3 = S4 = 1/4. Figure shows 
(red) the preimage under exp of x : [0, 1] — > S 3 with observations xq = (1, 0, 0, 0), 

x x = (0.1304,0.7923,0.4574,0.3821) 

x 2 = (0.5809,0.0381,0.3385,0.7393) 

x 3 = (0.5523,0.6251,0.5513,0.0172) 

x A = (0.2810,0.1241,0.6817,0.6640) 

made at to = 0, t\, ti, £3 and t&. Using Mathematica's FindMinimum on Y2k=l \\-BL,k\\ 2 > 
considered as a function of A, the best estimating geodesic from xo is found to have in- 
finitesimal generator Al = (1.40398,0.196766,1.05334), with preimage shown in black in 
Figure [3| We find 




Figure 3: The Conditional Minimum and Optimal Prior in Example 11 



B L> i = (2.7669,2.3129,2.0736) 

B Lt2 = -(1.0075,4.2867,1.4298) 

B Lt3 = -(2.1680,-3.6097,2.4150) 

B LA = (0.4086,-1.6359,1.7713) 

with sum in accordance with Corollary [6| When, as presently, observations tend to 
contradict the hypothesis that x is an integral curve of a field near A, the usefulness of 
the optimal interpolant is questionable. For less contradictory observations we find the 
preimage of x in E 3 more nearly piecewise-affine, consistent with Example [J]. □ 
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Corollary 7 Take n = 2 in Corollary^ If (x,A) is optimal then e si L is the midpoint of 
a geodesic joining y\ and y2- According as t G [to, or t £ [ti^tz], 

x(t) = e^* ^ 1 e (*-*o)^ x or 
x (t) = e fr-t)Bz.,i e {t-t x )A L xi 

where e SlBL ' x = y\e SlAh . 

Proof: By Corollary |4j for an optimal field A, e slBL > 1 e slAL = y\ and e~ SlBh ' 1 e slAL = y2 ■ 
Then s G [— s\, s%] h-> e sBL - 1 e SlAL G G is a geodesic from 2/2 to y\. □ 

7 Symmetric Priors on Sphere and Hyperboloid 

Define a nondegenerate symmetric bilinear form ( , ) on M 3 by 

{(vi,v 2 ,v 3 ),(wi,W2,w 3 )) = VlWl +v 2 w 2 + av 3 w 3 

where a = ±1. According as a is 1 or — 1, ( , ) is the Euclidean metric or the Lorentz 
metric. In either case ( , ) restricts to a Riemannian metric of constant curvature a on 

M a = {y GM 3 : (y,y) = a}. 

Then M + is the unit 2-sphere S 2 with the standard Riemannian metric, and M_ is the 
2-sheeted unit hyperboloid H 2 [7j. 

For vector fields X, Y on M a and a C 1 vector field A on M a , 

Oa,x(Y) = {V X A,Y)-{V Y A,X) = (X(A),Y)-(Y(A),X). 

Define vector fields B and C on M a by B(x) := (— x 2 ,xi,0) and C(x) := (0,0, 1) — x 3 x. 
Since C is the gradient of the height function, Qc,x = 0- Also 

B ,x(Y) = 2{X l Y2-X 2 Y l ) = 2{(-X 2 ,X 1 ,Q)-a(-x 1 X 2 + x 2 X 1 )x,Y). 

So 9g T x = 2((-A 2 , Ai,0) - a(-xiX 2 + x 2 Xi)x). Given C 1 functions (3, 7 : M CT -> R, set 

A := + p|| 2 = (/3 2 + 7 2 )(x? + x^). 

By(§,^ = /3^,x+^(/3) J B T -(S,X)ci/3 + A(7)C T -(C,A)d7. 

Let ^4 be invariant with respect to rotations of the first two coordinates of M a , namely 
j3 = (3 o pr 3 and 7 = 70 pr 3 where f3, 7 : R — > R and pr 3 : M CT — > R projects to the third 
coordinate. Then 

^grad(P|| 2 ) = ((^ 2 +7 2 )*!, W 2 +f)x2, (^ , +77)(l-^))-(^ 2 +7 2 +(^ , +77 )x 3 )(l-x 2 )x. 
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Also e~ A T x = (39z T x + X((3)B-(B,X)(dp)- T = 
2(3((-X 2 ,X 1 ,0) + wx) + i3'((-x 2 X 3 ,xiX 3 ,-w) +wx 3 x) 
where w := a (xiX 2 - x 2 Xx). Taking X = zW for a C 2 curve x : K — > M a , equation is 



equivalent to 

xf ] +a\\x^\\ 2 xi = OS 2 + 7 2 )x^i- 00' + 7f)(l-x|)x 3 xi + 

-2 ( 94 1) - /^W^ + ( 2 /3 + P r x 3 )wxi (6) 



x 



( ~ 2) +a||x( 1 )|| 2 x 2 = (,3 2 + 7 2 )x 2 ; x 2 - + 77) (l-xj)x 3 x 2 + 



2^ 1} + ^'xix^ + (2/3 + p , x 3 )wx 2 (7) 
4 2) +<t\\x {1) \\ 2 x 3 = -(^ 2 + 7 2 )(l-^)x 3 + (^/3 , + 77 , )(l-x 2 ) 2 + 

(2/3x3 - ~ xl))w. (8) 

By equation ([3]), for some 6 € K, 

allxWH 2 = (^ 2 + 7 2 )(l- a ; 2 ) + 6. (9) 

By rotational symmetry and Noether's theorem, 

w = o-(xi4 1} -x 2 x ( \ ) ) = j3(l-xl) + c. (10) 

where c G E. Set r(i) = \J x\ + x 2 = -y/ <r(l — x 2 ). For r|(oo, ai) nowhere-zero, choose 
■0 : (a0) a i) -> R so that xi = rcos^ and x 2 = rsin^. Then w = ar 2 ip^ and, by (10), 

^(t) = P + ^- 2 . (11) 
1 x 3 

So X3 : (ao,ai) — >■ K determines x|(ao,ai) up to quadrature. 
Proposition 1 fj||x (1) (t) - yl(x(t))|| 2 = b - 2c/3 + 27(1 - 2x 3 )x^ } . 

Proof: \\x^(t)-A(x(t))\\ 2 = - ^(-x 2 , x l5 0) - 7(-x 3 xi, -x 3 x 2 , 1 - x|)|| 2 = 
||x (1) || 2 + cr(/3 2 + 7 2 )(1 - x\) - 2ofiw + 27(x 3 xixj 1) + x 3 x 2 x { 2 l) + a(l - xjjx^) 

= cr6-2o-/3c + 27(-crx|x^ 1) +cr(l-xl)x3 1) ) = ab - 2acf3 + 2aj(l - 2x 2 3 )x ( 3 1] 
by (0), (10), and because x^xi + x^ 1 ^ + crx^ 1 ^ = 0. □ 
Corollary 8 Let /3 and 7 be constant. Then 

a J \\x w (t) - A(x(t))\\ 2 dt = (b - 20c)t + 2j(x 3 - 2x^/3) + c 

where cq is constant. □ 
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Proposition 2 Let (3 and 7 be constant. A solution x : R — )■ M a of ([?[), Q), 23 is 

constant if and only if, for some ipo 6 R, either 

1. x(t) = (cos(wi + ■i/'o)) sin(o;t + tpo), 0) /or a// 1, where U) 6 R is arbitrary, or 

2. 7 = and = (y/o-(l - h 2 )cos((3t + ifj ), ^/a(l - h 2 ) sin(/3t + tp ),h) for any h 
with cr(l - /i 2 ) > 0. 

Proof: Suppose x^(t) = h for all t. If h 2 = 1 then x is constant at ±(0,0, 1). Conversely 
these constant solutions satisfy (p|, ff\, ph. Also 6 = c = 0. Suppose now /i 2 / 1 and 



/3/0. By (11) 



x(i) = (y/ a(l — h 2 ) cos(wt + tfjo), \J cr(l — h 2 ) sin(wi ± ^q)) h) 
- C I 1 I 1 

where u) := (3 H — r and ^0 £ R- Equations (6), (7) are together equivalent to 

1 — /r 

c 2 + 2(l-/i 2 ) 2 /3c-((l-2/i 2 )7 2 + 6)(l-/i 2 ) 2 = 0, 

and equation |8j is equivalent to (27 2 (1 — h 2 ) — (2(3c — b))h = 0. So either h = and 
( c + 6)2 = ^2 + ^2 + b) or 2/ 3 c = 27 2 (1 -h 2 ) + b and c 2 + 7 2 (1 - h 2 ) 2 = 0. In the first case, 
c is arbitrary, depending on b. In the second case c = = 7, and h with <r(l — /i 2 ) > is 
arbitrary. □ 

Example 12 Lei /3 and 7 &e constant. If /3 = = 7 £/ien a// constant curves x satisfy 
@; If P ¥^ ® an d 7 = £/ien £/ie on/y constant solutions are x(t) = ±(0,0,1) 

and x(t) = (cos V'O) sin-^O) 0) where ipo E R. If neither j3 = nor 7 = 0, £/ie on/y constant 
solutions are x(t) = ±(0,0, 1). □ 



Substituting for ||xW|| 2 and w in (j8|), we find 

x ^ + bx 3 = - 27 2 (l-x 2 )^3 + 77 / (l-^) 2 -/S / (l-x 2 )c±2^ 3 c 
and, integrating with respect to X3, for some d G R, 

(4 1 ') 2 + ^3 = 7 2 (1 - xl) 2 - 2^(1 - x|) + d. (12) 
More can be said for particular choices of j3 and 7. 

7.1 7 = and /3 Constant 



Let x : R — > M a be a solution of (|6j), (|7j), (|8|), where 7 = and j3 is constant. By Proposition 
[TJ the integrand for J(x, A) is constant. If X3 is constant then x is given by Proposition [2] 
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Theorem 6 Let x 3 be nonconstant. For some A with a(l — A 2 ) > and some Vq G E, 
according as a = ±1, 

x(i) = \Jl - A 2 sin 2 (te + t;o) (cos^(i), sin ip(t), 0) + A(0, 0, sin(te + « )) or 

x(i) = ^/A 2 cosh 2 (ie + v ) - 1 (cos^(i), sin^(i), 0) + A(0, 0, cosh(ie + u )) 
where, for some o\ = ±1, either A 2 = 1 and ip(t) = t(3 + ^o? or cr(l — A 2 ) > and 

ip(t) = Vo + tp db arctan ^a/ 1 — A 2 tan(ie + i>o)^ + 2to7t where m G Z, or 

^(i) = + *0 ± arctan f ^^= , 
according as a = ±1. Conversely, any such x satisfies |6p, ([?[), |S|). 



Proof: By Q, (x 3 1} ) 2 = (20c - 6)s§ - 20c + d. If 6 = 20c, then (x 3 1} ) 2 = d - 20c, and 
(i = 20c = 6 since otherwise x 3 maps onto the whole of E. So x 3 is constant when b = 20c. 

For b / 20c, we claim that, for some A with a(l - A 2 ) > 0, (x^) 2 = (6 - 20c) (A 2 - x 2 ) . 
This holds trivially if (20c - b)x 3 (t) 2 - 20c + d = for some tel. If (20c - 6)x 3 (t) 2 - 
20c + d ^ for all t, then x 3 is strictly monotonic, and x 2 is bounded above or below 
according as a = ±1. So lim^-too x 3 (£) exist, lim^-too X3 (i) = 0, and 

= lim x^\t) 2 = (20c- b)\l - 2(3c + d where A± := lim x 3 (t) 

t— >±oo " t— >±oo 



which proves the claim. Set e := y\b — 20c| > 0. Then (x^) 2 = (b — 2c0)(A 2 — x§). 

• Suppose 6 > 20c. Then x 3 (t) = ±Asin(te + Vo) where vq G E. If a = —1 then 
£3 (t) 2 > 1 for all t. So a = 1 and A 2 < 1. Then 

x{t) = \Jl-\ 2 sin 2 (re + f ) (cosV>(*), sin^(t), 0) ± A(0, 0, sin(ie + v )). 



Differentiating, and substituting for ijj^ from equation (11), 



,.,.1)112 _ e 2 A 4 sin 2 (te + -i> )cos 2 (te + ?;o) | r 2 v 2,,, 
1 - A 2 sin 2 (ie + v 



+ e z X z cos z (te + vo) + 



(1 - A 2 ) sin 2 (t e + v )(P + - , 2 . 2 J 2 . 

1 — X z sin (ie + vq) 



Equating this with the right hand side of ([9]), c 2 = e 2 (l — A 2 ). Then the expressions 
for ip follow on integration of (11). Substituting for x in ^ , Q , Q , these equations 
are indeed satisfied. 
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□ 



Suppose If b < 2/3c. Then, for some vq £ R, x$(t) = ±Acosh(ie + vq). If a = 1 then 
xslt) 2 < 1 for all t. So a = -1, A 2 > 1, and 



x(t) = yA 2 cosh 2 (te + w ) - 1 (cosV'(t),sin^(t),0)±A(0 ) 0,cosh(ie + 'Uo))- 
Proceeding as before, c 2 = e 2 (A 2 - 1), and Q, (Ft), © are verified on integration of 



(11). 



Example 13 Figure [^] shows (red) x : [0,1] — > S 2 satisfying |fi]), |?|) ; a = 1, 

P = —1, 7 = ; A = 1, e = 7, vq = ipo = 0. Figure^ also shows (black) the orbits of A, 
and (blue) a numerically calculated conditional minimum x a : [0,1] -> S 2 , toii/i 

^o(O) = x(0) = (1,0,0) and x Q (l) = x(l) = (cos 1 cos 7, — cos 7 sin 1, sin 7). 

Then J(x a ,A) ~ 0.519 and J(x,A) = 49. Repeating with e = \/2 gives x ~ x, J(x D ,^4) ~ 
2.011 and J(x) = 2. 




Figure 4: A Non-Minimal Conditional Extremum (red), with the Minimum (blue), in 



Example 13 



Figure^ shows x : [0, 10] — > S 2 with a = 1, (3 = -I, 7 = 0, A = 1. 

Figure^ shows x : [0, 10] — > S 2 with a = I, (3 = —1, 7 = 0, A = 0.5, e = 1.1, vq = ipo = 0. 
□ 

Example 14 The Poincare disc D is the open unit disc in IR 2 with Riemannian metric 

{v\wi + V2W2) 



(v,w) := 4- 



l-z 2 + z 2 ) 2 
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Figure 5: a = 1, (3 = — 1, A = 1, e 2 = 1.1, v o = tpo = in Example 
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where v,w are tangent to D at z. Geodesies of D are either circular arcs orthogonal to 
the unit circle bounding D, or diameters. An isometry to D from the path component of 
(0, 0, 1) in H 2 is given by 

v s (2/1,2/2) 

y ^ <J)(y) ■= . 

1 + 2/3 

Rather than directly plot critical curves in the hyperboloid H 2 , a clearer picture is obtained 
by mapping to D. Figure^ shows (f)ox : [0, 10] — > D (red) with a = — 1, f3 = — 1, A = 1.01, 
e = 4, vo = and ipo = —ir/3. Integral curves of the field corresponding to A are shown 
as circles (black). Where the conformal factor in the metric is large, towards the end, the 
conditional minimum x tends to agree with the orbits of the prior field. Disagreements with 
the prior field are greater towards the beginning. □ 




Figure 7: a = —1, /3 = —1, A = 1.01, e = 4, vq = 0, ipo = —ir/3 in Example 14 



7.2 f3, 7 Constant and 7 ^ 



Let x : R — > M a be a solution of (|6j) , ([7|), ([8|), where /3 and 7 are constant with 7 7^ 0. By 
Theorem [3], A is reflexive if and only if j3 = 0. If X3 is constant then x is given by part 1 
of Proposition [2] 

Theorem 7 // X3 is nonconstant then, for some a E C 



x 3{t) = ±yp(~ft + a;g 2 ,g3)-o~ 

where 5 := (2/3c— b— 2j 2 ) / (3*/ 2 ) , p is the Weierstrass elliptic function and, for some del, 

52 = l25 2 + d 
93 = -8P-M. 
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Proof: By H, (x^) 2 



f(x\ + 35x1 - d / 4 ) where d ■= ~ 4 (7 2 - 2 / ?c + d )h 2 - Setting 



y(s) := x 3 (s/7) + 6, we find ( — j = 4y -52^-53. □ 



So (11) gives x up to quadrature in terms of a known periodic function X3 : 



Example 15 Take a = 1, /3 = and 7 = 1. Since A is conservative it is reflexive, by 
Theorem [3J and i/ie reverse x of a conditional extremum is also a conditional extremum 
for the reversed data. Figure^ shows (red) the solution x : [0,14] — > S 2 of |fi]) ; |?|) ; |#|) 
with x(0) = xq = (\/3/2, 0, 1/2) and x^^O) = (0,1,0). TTie integral curves of A, namely 
the gradient of the height function, are shown in black. We find a ~ 1.14811 + 1.74899?, 




Figure 8: A Conditional Extremum for the Conservative Prior Field in Example [l5| 
b = 0.25, c = -0.5, g 2 = 4.75, g 3 = 1.875. □ 

Example 16 Take a = — 1, (3 = —1 and 7 = 2. Figure^ shows (red) (p o x : [0, 1] — > D 
where x : [0,1] — > H 2 is the conditional extremum with x(0) = xq = (0.1, 0.1, Vl-02), 
x(l) = xi = (8.99009, -7.34992, 11.6552) and x (l \0) = (0.9^/1-002,0,0.9). Integral curves 
of the prior field A spiral inwards (black). Numerical calculations suggest x is a conditional 
minimum. We find g 2 = 1.17393, g-$ = —0.220814 and a ~ —2.08599. The prior field A is 
effectively contradicted by the supposed observations xq and x\. The conditional minimum 
x makes the best of difficult circumstances. 

On the other hand, if the order of xo and x\ is interchanged, these might well lie on 
an integral curve of a vector field near A. Figure [5] shows the corresponding conditional 
minimum (blue, dotted). Disagreements with A tend to be concentrated near the centre 
of the disc, where the conformal factor in the Riemannian metric is relatively small. By 
Theorem^ A is not reflexive, because f3 7^ 0. □ 



Acknowledgement: This paper is dedicated, with great respect, to the memory of Pro- 
fessor Jerrold E. Marsden. 
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Figure 9: Conditional Minima for the NonConservative Prior Field in Example 
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